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I. INTRODUCTION 

The relevance of studying the Kondo lattice model 
(KLM) has not decreased. It is believed to be at the heart 
of the physics of both the heavy fepmion materials, in 
its antiferromagnetic (AFM) versionJH and also the man- 
ganites, when local moments and conduction electrons 
interact via the ferromagnetic (FM) Hund's coupling.El 
In the first case, the well understood behavior of the sin- 
gle (or few) impurity Kondo model permeates much of 
our current understanding. However, the interplay be- 
tween the local Kondo physics and the non-local RKKY 
interaction in a lattice environments remains elusive in 
current approximate schemes, although it mav nif|'li a y a 
prominent role close to quantum critical pointsQO'Q'LJ or 
even away from them. In this respect, a more thorough 
understanding of the one-dimensional (ID) case might 
be fruitful, even in light of the peculiarities of ID sys- 
tems. Furthermore, the study of the ID KLM is impor- 
tant in its own right for the analysis of some quasi-one- 
dimensional |-oiganic compounds such a& (Per)2M(mnt)2 
(M=Pt,Pd) ! Ma and (DMET) 2 FeBr 4 .B 

A fairly complete ground *t&tje phase diagram has been 
established for the ID KLMBll 2 ] On the one hand, the an- 
tiferromag-netic model at half-filling has both charge and 
spin gapsilia For lower band fillings there is a quantum 
phase transition frpm a paramagnetic ground state to a 
ferromagnetic oneJl3 On the other hand, the ferromag- 
netic model at half-filling is also insulating with a Hal- 
dane type spin gap£j For lower band fillings, the numer- 
ical evidence shows three distinct phases: a phase with 
partial ferromagnetic order and incommensurate spin 
correlations, a fully saturated ferromagnetically ordered 
phase and a region with phase separation where two kinds 
of ground state seem to compete. The energy scales 
where the transitions take place for both models are of 
the order of the Fermi energy. Finally, there is strong nu- 
merical evidence in favor of a Luttinger Liquid behavior 
in the paramagnetic phase of the AFMirKLM, even for 
considerably large coupling constantsJlail^'U Such phe- 
nomenoloawjis beyond a simple RKKY versus Kondo type 
of picturejla as proposed by Doniach for the higher di- 



mensional models .□ In fact, the level crossing responsible 
for the quantum critical point (QCP) is related to long 
wavelength modes. This is in contrast with the short 
wavelength spin modes involved in the paramagnetic- 
antiferromagnetic transition that the-JDoniach's picture 
envisages. As pointed out beforejl^lla the missing ele- 
ment is the lowering of the conduction electron kinetic 
energy with the alignment of the localized spin as in the 
double exchange mechanism. Essentially, this is the rea- 
son why the ID FM and AFM KLM have similar phase 
diagrams. 

Despite these successes, it would be considerably 
more informative if some analytical understanding could 
be gained. .Even J-hough there have been some par- 
tial successesliaocJ'EJ, there is still room for improve- 
ment. Motivated by the enormous success of renor- 
malizati(j»|^KQTOr-(pG) analyses in the few impurity 
problem jHdO'ESE^Ej we set out to apply scaling ideas to 
the ID lattice case as well. However, an RG treatment 
of the KLM has never, to our knowledge, been achieved. 
Technically, although we know how to progressively dec- 
imate the spins or the conduction electron states, no one 
has devised a way of doing both simultaneously, specially 
with the incorporation of local Kondo physics. 

It is the aim of this article to put forth such a deci- 
mation scheme for one-dimensional models of spins and 
fermions, in particular the anisotropic Kondo lattice 
model. We draw a great deal of inspiration from the orig- 
inal work of Anderson, Yuval and Hamman for the single 
impurity Kondo modeljHaO'Efl mapping the KLM into a 
classical Coulomb gas, which is then decimated by stan- 
dard methods Eli This task is made simpler by the use of 
bosonization methods. We therefore study the stability 
of the non-interacting ground state with respect to the 
Kondo interaction as a function of the coupling constants. 
Our study reveals that there is no "weak coupling" flow 
in the entire parameter space. Nevertheless, the differ- 
ent "strong coupling" flows of the RG equations allow us 
to assign the magnetic ground states that emerge, estab- 
lishing the phase diagram for both signs of the coupling 
constant in a unified fashion. W,hil<&i|e*p| approach in 
part builds upon previous studies ^tBESl 21 ! it also puts 
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on a firmer basis the procedure of neglecting backward 
scattering terms in the KLM away from half-filling. As 
another application of our Coulomb gas treatment, we 
also establish the phase diagram of the one dimensional 
Ising-Kondo model.c3 

In Section II, we develop a path integral formulation of 
the bosonized ID KLM. The partition function is mapped 
into a two-dimensional generalized classical Coulomb gas 
in Section III, In Section IV, the RG equations of the 



Coulomb gas are derived and solved. Their physical in- 
terpretation is given in Section [v], where an effective 
Hamiltonian for the renormalized Coulomb gas is ob- 
tained. The phase diagram of the model is established 
in Secti on (y| . The Ising-Kondo model is discussed in 
Section [VII| , where its phase diagram is established. We 
wrap up with a brief discussion of th e rel ation between 
our and previous results in Section VII] . Some more 
technical developments can be found in the Appendices. 



II. PARTITION FUNCTION 

We start by writing the ID KLM Hamiltonian. The 
traditional Kondo model is isotropic in spin space. Since 
we are going to use Abelian bosonization, it is natural to 
break the SU (2) symmetry down to U (1) 



a) 



where Cj a destroys a conduction electron in site j with 
spin projection a, Sj is a localized spin \ operator and 

Sj = lE^c^crc^, the conduction electron spin den- 
sity. We will focus on the continuum, long-distance 
limit of the conduction electrons. In this case, one 



J 



can linearize the dispersion around the non-interacting 
(J z ,± = 0) Fermi points ±kp, where kpa — and n 
is the conduction electron number density, and take the 
continuous limit of the fermionic operators in terms of 
left and right moving field operators^ 



H 



-ivf 



E / dx {A 



a, J : 



3 a,0,s 



i a,/3,s,s' 

where a, /3 = L or R, vp = 2tsmkpa is the Fermi ve- 
locity and a is the lattice spacing. The field operators 
can now be bosonized with theJnclusion of the so-called 
Klein factors in usual notationca 



\l2na 

_ F L : a p -i^[4 >cr (x)+e tr (x)]-ik f 



One can then rewrite the Hamiltonian in terms of the 
charge and spin fields 



as 



with: 



e ,.(x) = (6 1 (x)±e l (x))/V2, 



H = H Q + Hi +H{ + H b z + H b ± , 



(2) 



Ho = y E f dx(d x M 2 + (d x e„y 



H 



H f z = E Jl\J-9«d x <p s {x)S z {x) 

X 

y -ieV ^ 8s(x) cos \j2^<j) s (x)} S- (x) + h.c. 
z — ' 2na I 

X 

2 jb r 

H b = } — - sin \J2Trg p (f) c (x) + 2k F x] sin J2%g a (j> s (x) S z (x) 
z — ' net L 

X 

b 



Hi = 



E 2~ e ^ ^ cos [V^fp^c (x) + 2k F x] S- (x) + h.c. 



(3a) 
(3b) 
(3c) 
(3d) 
(3e) 



where g a 



1 and a ~ k F . Hq is the free 



2 . We have introduced the new parameters g a and 



bosonic Hamiltonian written as a function of 6 S c and 
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g p for future use. The "relativistic" description enforced 
by us broke the interaction term in two different compo- 
nents: forward-scattering, H* , and back-scattering, H b . 
They involve the spin current and the 2k F component 
of the magnetization of the non-interacting electron gas, 
respectively,E3Eil since it is well known that the main con- 
tributions to the spin susceptibility of the electron gas at 
low frequencies come from q ~ and q ~ 2k F . For 
further generalizatiomwe will consider the 4 parameters 
J±\, as independento It is important to note that the 
cosines and sines of the bosonic fields in Eqs. (||) are just 
a short form notation. Forward and backward Klein fac- 
tors do not have common eigenvectors.^ Thus, we shall 
not neglect their contribution to the simultaneous treat- 
ment of H* and H b . 

A quantum system of dimension d cap-be mapped into 
a classical system of dimension d+ The single im- 

purity Kondo problem has effective dimensiapcf-^p 0. The 
work of Anderson, Yuval, and HammanEd'EBEil showed 
that it can be mapped into a d = 1 classical Coulwnbngas, 
where the extra dimension is the imaginary time.EMa We 
will extend this idea and map the ID KLM into a d = 2 
classical problem. As usual, the starting point is the par- 
tition function: 



Z = Tr 



o -p(H +H I +H b ) 



(4) 



We will rescale the Hamiltonian and (3 by the Fermi veloc- 
ity. This introduces the dimensionless coupling constants 
Jz,± = aJ J F ± as we ll as /3 = vp(3. Following the standard 
prescription we divide (3 in infinitesimal parts, 



Tr 



-8tH 



In order to proceed to a path integral formulation we 
choose the S z basis for the local moments and the coher- 



ent states for the bosonic fields 00 The next step is to 
introduce an identity resolution between each exponen- 
tial in the product 



where we used |C) to denote a general vector in the basis. 
We now expand the exponentials in powers of St, 



Z 



n=0 



(ST) n {<;(j)\H n \UT j+ i)) 



There are two possible spin configurations for a given 
pair of consecutive instants along the imaginary time di- 
rection: 

1. if there is no spin flip between them, the only con- 
tributing terms are from the z components of the 
Hamiltonian (3b) and (3d) 



(s(x, t + St)\ S z (x) d x cj> s (x) \s(x, t)) ; 

(s(x, t + St)\ S z (x) sin [y/2irg p <f) c (x) + 2k F x] 



sm 



\]2ng a 4> s (x) \s(x,t)) 



(5) 



2. if there is a spin flip between these 2 instants, then 
only (pScj) and (|33) contribute 



(s(x.t + Sr)\e V <>° * w , 
S~ (x) +h.c. \s(x,t)) ; 



(s(x, t + 5t)\ e 



cos 



y/2irg cr (t) s (x) 



\^/2Tig p 4> c (x) + 2k F x] 



S (x) +h.c. \s(x,t)) . 



(6) 



The two possible processes above are illustrated in Fig. 0. 
Several bosonic operators can fit inside the example given 
in that figure, for example 



4r ( \Z5" Pz I F h F ^ F k F ^ sz ( x > T 3) s+ («. T2) S z (x, n ) 



d x 4> s (x,ti) 



r 



where we regrouped terms to separate Klein factors, local 
spin operators and bosonic fields. 

The lattice parameter in the Euclidean-time direction 
is set by the bosonic cut-off: St = 2na. Keeping the 
leading order terms, we write the partition function as 
the sum over all d = 2 Ising spin configurations of the lo- 
calized spins, Klein factors and a functional integral over 



the bosonic variables. Since there are different types of 
bosonic exponentials ("vertex operators"), coming from 
the different interaction terms in Eq. (||), we now intro- 
duce new Ising variables which we will call "charges" in 
order to do the bookkeeping. They give the sign of the 
corresponding bosonic field in the accompanying expo- 
nential according to the following scheme: 
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Figure 1: Example of a spin history. 1 stands for S z operators 
(Eq. (|)) and 2 for Stones (Eq. (|)). 



fugacity 


(m,e,c) 


Number of particles 


yi = Jin 


(±1,±1,0) 


Ni 


V2 = Jl/2 


(±1,0, ±1) 


N 2 


yi = Jin 


(0,±1,±1) 


N 3 



Table I: Particles in the ID KLM and their charges. 



1. m (x, t) — S z (x, t + 5t) — S z (x, r) = ±1 gives the 
sign of 6 S (x, t) (Eqs. @ and @) 



2. e (x, t) = ±1 gives the sign of 4> s (s,t) (Eqs. ( |3c| ) 
and @) 



3. c (x, t) — ±1 gives the sign of C (a;, r) (Eqs. ( |3d| ) 
and (pq)) 



Note that only m (x,t) is always tied to a localized spin 
flip process, its value giving both the sign of the 9 S co- 
efficient and the change in S z . With this notation, each 
point in the Euclidean "space-time" is labeled by a triad 
of values (m, e,c). We call a "particle" a point where 
(m, e, c) 7^ (0, 0, 0). Each kind of particle matches a cer- 
tain incident in the history of a spin. From Eqs. ||, we 
can read off the existence of three breeds of particles, each 
one with its respective fugacity. Table | summarizes the 
notation that we will use. 



Denoting r/j as the space-time position of particle j of 

2, J 

the partition function as: 



type I = {1, 2, 3} and T>r\ = J\i=i Yifli ^Vj > we can write 



E E E v^ve^v^ 

r _i at. at. at. n r _ _i ** 1 



IVi N 2 N 3 

Vi Vi Vz 



{a} N 1 ,N 2 ,N 3 =0 {m,e,c} 



\N 2 \N 3 \ 



(Klein factors) 



exp i -S - J{J-ga ^2 j drd x <p s (x, r) S z (x, r) + iJ — ^ ^ 

[ x » 5cr |=1,2 ,j! 



JV 3 



E E e ^ Wi) + E E l c (4) & + 2 ^ c W) 4] 



(7) 



'=1,3 ^ 



'=2,3 ^ 



where {c} stands for all Ising spin configurations, 
{m, e, c} represents all possible sets of J^i^i particles, 
Xj is the space coordinate of particle r) l j and Sq is the 
free_Gaussian bosonic action in the variables 4> s ^ c and 
6* s , c EiJ. There are several restrictions over {m, e, c|,-|USu- 
ally called neuteality conditions in the bosonizatiorO and 
Coulomb GasEZl literatures. In the ID KLM they are 
more stringent than usual, ensuring the compatibility of 
the sums over {a} and {m, e, c} in Eq ([?]). Therefore, we 
will call them strong neutrality conditions (see Appendix 
^ for the derivation of these conditions). They are: 

1. for each space coordinate the m charges must be 
neutral, Y,i m i x 'fixed, n) = 0, 

2. the total e charge must be neutral, J2i e ( x i> T i) = 0) 

3. for each space coordinate the total charge c must 
be an even integer J2i c ( x fixed, T i) = 2n, n e Z; 
and the total charge in the entire space-time must 
be zero, J2i c ( x i, T i) = 0. 



We can immediately see two consequences of these con- 
ditions. The most obvious is that the sign of jf' 6 is 
irrelevant since, from condition 1, the total number of 
spin flips in the time direction N\ + N2 is even. The 
other consequence is more subtle and more surprising: 
the complete cancellation of the Klein factors and the 
product of a z (7?!) in Eq. (0), 

N 3 

(Klein factors) J\ <? z (Vj) = 1- 

j=i 

This result plays a central role in the renormalization 
group treatment-,of a single Kondo impurity in a LL by 
Lee and Toner .E3 Moreover, it leads to: 



2k F c (*) Xi = < , (8) 

i I Ak F aI, I € Z 
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in each contribution to the partition function. The 2kp 
terms appear whenever there are particles of type 2 and 
3 (see the definition of the charge c and Eqs. (3d) and 
(|3e|))- Due to their oscillatory nature, configurations with 
these particles will be strongly suppressed in the statis- 
tical sum and the corresponding terms (with fugacities 
j/2 = J±/2 and ?/3 = Jz/2) will be irrelevant in the RG 
sense. This irrelevance criterion is precisely the same as 
the one used from neglecting Umklapp scattering awajz. 
from half-filling in models like the Hubbard model.Ej 
However, we stress that the situation here is far less 
trivial than in the Hubbard model, since we have both 
fermions and spins, and the latter have no independent 
dynamics, thus hindering a rigorous analysis (an inmar- 
tant exception to this is the Heisenberg-Kondo modelKa). 
In our treatment, on the other hand, spins and fermions 
are treated on the same footing and lose their indepen- 
dent identity. After the mapping to a Coulomb gas, the 
irrelevance criterion becomes identical to other models 
where its applicability is firmly based. We have thus es- 
tablished a more rigorous basis for neglecting the back- 
ward scattering terms in the ID KLM, as has been done 
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by Zachar, Kivelson and EmeryEi 

This situation changes when the conduction band is 
at half filling. In this case, Akpa = 27r and these terms 
disappear from the effective action, making all particles 
equally probable. This commensurability condition is 
similaE-tp the one for the Umklapp term in the Hubbard 
model.E3 It is interesting to note that only the combina- 
tion Akpa appears in our formulation. Since we bosonized 
the non-interacting conduction electron sea, we must use 
kpa = irn/2, leading to 4fci?a = 2irn. Even if foip-some 
reason a large Fermi surface should be considered, e£l this 
would not change any of our results since for a large Fermi 
surface Ak F a — 2ix (n + 1) = 4fc^amod27r. 



III. COULOMB GAS 

The bosonic fields in Eq. (0) can now be integrated out, 
partially summing the partition function. The result can 
be understood as an effective action for the spins and the 
variables (m, e, c) 



S, 




E 



dndr 2 C -?^S*(l)S*(2) 

r 12 

exp (ie (n) (p) 



V E E / dr ^^W' S* (x, r) 

n x J 

E ~~f^ ( m (™) + e ( n )) ( m (p) + e (p)) + ( m ( n ) ~ e ( n )) ( m (p) _ e (p)) 

n>p 

In | r np | c(n) c(p) + 2ikp ^""^ c (ra) xi + short range interactions, 

n 

I 



(9) 



where 

Tjk = Tfe - Tj , 

Zjk = Xjk + iTjk = r jk e lVjh . (10) 

In addition to the long range universal interactions, 
this procedure also gWes rise to short-ranged terms that 
are cutoff dependent Jiil H}ese are similar to those found 
by Honner and Gulacsili3E3 by smoothing the bosonic 
commutation relations. In contrast to their treatment, 
though, here they are a manifestation of the bosonic field 
dynamics. Following Zachar et. aljtB we will focus on the 
universal long range part of the action and neglect these 
terms. 

Upon integrating by parts in imaginary time, spin time 
derivatives become the charges we denote by m. Finally, 
we can rewrite all long range terms in the form of a gee-, 
eralized CG action in two-dimensional Euclidean spacecll 



as 

N 1 ,N 2 ,N 3 =0{m.e.c} J 

(11) 

with 

S eff = E 1 _ ln I m 171 fa ) 

+ 9a In |r y I e (77,) e (r)j) + g p ln |r y - 1 c (%) c (%) 

- inipij [e (rji) m (r]j) + m (rji) e (rjj)] | 

+ 2ik F ^2c(rn)xi, (12) 

i 

where « = 1 — J//tt. In the above effective action we 
have dropped the superindex indicating the particle type 
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in order to unclutter the notation. It is now unneces- 
sary as the dependence with the history of a spin has 
disappeared. 

In most other similar CG mappings, the coefficient of 
the term in 4»™ is an integer and goes by the name of con- 
formal spin£il Then, the ambiguity of 2irl, I G Z in the 
angle is irrelevant. In this case, however, k can assume 
non-integer values. What guarantees that the theory is 
actually well defined is the strong neutrality condition 
1, which leads to a cancellation of the Riemann surface 
index /. 

The integration by parts that we performed is equiv- 
alent to applying the duality relation d x (f> s — id T 9 s to 
Eq. (0), integrating by parts and then tracing the bosonic 
fields. Alternatively, at the Hamiltanian level, it is also 
equivalent to applying the rotation!^ 
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RG step 
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Figure 2: Length rescaling in the CG. Originally distinct 
charge configurations are identified at the new scale. 



u 



(13) 



to Eq (g) before going to a path integral and tracing out 
the bosons. Hence, there is a strong link between ows-CGi 
formulation and previous results in the literature jlallSEil 
More importantly, the interpretation of our results should 
be understood in this rotated basis that mixes spins and 
bosons. 

The effective action in Eq. ( |l~2] ) can be viewed as de- 
scribing the electrostatic and magnetostatic energy of 
singly charged particles with both electric and magnetic 
monopoles. These satisfy electric-magnetic duality in the 
sense that the action is invariant under the exchange 

2 

e ^ m and g a ^ while k is unchanged. This is 
analogous to the Dirac relation between electric and mag- 
netic monopoles. Furthermore, these particles possess a 
third electric- like charge (c), unrelated to the two pre- 
vious ones. The partition function sum now has been 
reduced to considering all particle configurations, blend- 
ing spins and bosons in this Coulomb gas representation, 
where we have particles plus neutrality conditions. 

A partially traced partition function allows us to link 
problems that are originally quite distinct. For exam- 
ple, the only difference between the CG's of the single 
impurity Kondo problem and the problem of tunneling 
though an impurity in a Luttinger liquidcHl are the neu- 
trality conditions. Analogously, the two-channel Kondo 
problemc3'cil and the double barrier tunnelingpa can be 
mapped into each other with the same neutrality condi- 
tions. The KLM also has an unsuspicious counterpart 
in the li 
liquids L. 

is analogous to a spin flip process that scatters a boson 
from an up spin band to a down one and vice-versa. In 
particular, the two problems give the same effective ac- 
tion (with different neutrality conditions) if we disregard 
the backward-scattering terms in Eq. (g) and consider 
the anisotropic case k = 1 (<// = 0). 

In the following section, we will derive the Coulomb 
gas renormalization group equations following closely 



ire: two weakly coupled spinless Luttinger 
!£l The tunneling from one LL to the other 



the review by Nienhuis.0 As expected, the procedure 
strongly resembles the renormalizj^i©ji--gm««— aealysis 
of the tunneling between 2 LL'sMsoSMaEa The 
Coulomb couplings g a . p are equal to 1 for non-interacting 
conduction electrons. However, the same RG equations 
will apply to the case of conduction electrons with an 
SU(2) non-invariant forward scattering interaction. In 
this case, the initial values of-a^p are the corresponding 
Luttinger liquid parameters rB We will not dwell upon 
this case here, but its phase diagram is analogous to the 
one we will derive below. 



IV. RENORMALIZATION GROUP EQUATIONS 

The philosophy of the renormalization group is to 
sum the partition function by infinitesimal steps and 
find recursive equations for the coupling constants while 
keeping the same form of the effective action. In a 
Coulomb gas each step corresponds to three distinct pro- 
cedures: length rescaling, particle fusion and particle 
annihilation.^!! In order to implement these procedures 
all the particle fugacities must be small and we are forced 
to impose J h z and j{' b <C t. 

The first step consists of integrating large wavelength 
modes and then rescaling-parameters so as to reconstruct 
the original action form.E2l This corresponds to the overall 
length rescaling: 



l-dl 



(14) 



in the action and the partition function measure Dr\. 
Upon rescaling we lose the ability to distinguish cer- 
tain previously distinct charge configurations as exempli- 
fied in Fig. ||. Applying Eq. ( [l4] ) to the effective action 
Eq. ( |l^ ) and expanding the logarithm for di <C 1, we 
obtain: 
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fugacity 


(m,e,c) 


Number of particles 


G = 


(±2,0,0) 


N 4 


G = 


(0,±2,0) 


N 5 


r = 


(0,0, ±2) 


N 6 



Table II: New particles created upon rescaling and their 
charges. 



This equation gives the dimension of the corresponding 
operator and leads to the standard relevance criteria for 
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RG step 



S ef f - S, 



eff 



2 k 



— m (r)i) m (ijj) 
9c 



+ 9*e (Vi) e + g p c (rji) c (rjj)} dl. 

The neutrality conditions can be used to rewrite the last 
term as a single sum over sites: 



S e ff - S, 



eff 



— m (rji f 

9a 



+ 9ae(r]i) 2 +g P c(r]i) 2 



dl. 



(15) 



The integral over r] is the sum over all possible particle 
positions, and its rescaling leads to: 



dm 



(1 - d£)° 



(16) 



We have left the dimension d unspecified for the fol- 
lowing reason. Since a particle can only exist at the space 
coordinate where a spin exists, we can define two impor- 
tant limits in the KLM. If the Kondo spins are separated 
by a distance greater than dl, the sum over identical 
configurations is one dimensional (d = 1), as in the sin- 
gle impurity case (see (6) and (c) of figure ||). This is 
the dilute limit or "incoherent regime" of the Kondo lat- 
tice, where the scaling proceeds exactly as in the single 
impurity Kondo problem in a LL as found by Lee and 
Toner £3 In contrast, when dl is larger than the-distance 
between Kondo spins we are in the dense limhx£l or "co- 
herent regime" of the Kondo lattice. In the latter case, 
the identification of initially distinct configurations can 
involve charges at different space coordinates, implying 
that d = 2 (see (a) in figure ||). We will focus on this 
coherent regime and set d — 2 from now on. 

Collecting Eqs. @ and @ we can express the parti- 
tion function once again as a Coulomb gas by redefining 
the particle fugacities. A particle with charges (m, e,c) 
has its fugacity Ym,,e,c renormalized as 



dY„. 



Y 

1 n 



(17) 
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Figure 3: "Close pair" processes in the RG step: (a) particle 
annihilation and (b) particle fusion. 



bosonic operators. However, we have so far disregarded 
two possibilities. Suppose that a pair of initially distinct 
particles is within range, of the new smallest scale. Fol- 
lowing Anderson et alxA we call it a "close pair". After 
the RG step we can no longer resolve these two particles 
as separate entities. On the one hand, if the particles 
have precisely opposite charges we have a "pair annihi- 
lation" (see (a) in Fig (||)). The residual dipole polar- 
ization of this pair renormalizes the interaction among 
the other particles, leading to the RG equations for g a 
and g p . Note that n is an RG invariant. On the other 
hand, if the pair is not neutral the particles are fused 
into a new particle carrying the net charge (see (b) in 
Fig (||)). This last process may actually create parti- 
cles previously absent in the gas. There are three new 
kinds of particles created upon fusion in the dense limit 
with initial conditions g a ~ g_p ~ 1. Their charges and 
fugacities are listed in Tab. |fi|. These new entities cor- 
respond to originally marginal operators that are absent 
in the bare problem (their physical meaning will be dis- 
cussed in the next section). Other particles with higher 
charges could also be considered, but from Eq. (|l7| ) it is 
clear that they are highly irrelevant and therefore can be 
neglected. Collecting the annihilation and fusion terms, 
derived in Appendix [^, and adding the dimensionality 
equation (|17|), we complete the renormalization group 
equations. Away from half-filling, where the backward- 
scattering terms are irrelevant, particles with fugacities 
?/2,3 and T can be disregarded. Thus, only configurations 
involving the fugacities yi, G and G need to be consid- 
ered. On the other hand, at half-filling all particles from 
Tabs. | and [H] should be included. This leads to the fol- 
lowing renormalization group equations 
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Away from half-filling 



At half-filling 



dyi /„ 1 /k 2 \\ sin(27TK) 

2 - n — + 9* ))vi + s -2/1 G + G , 



d£ \ 2 \ j„ yy 2k 

rlC 

— = 2(l-g a )G + nyl 

* =2 ( 1 -^J G + ^' 

1 dingo- sin(27TK) / k 2 \ , k 2 ~ , ^ 9 
' --fl* + — G 2 -g CT G 2 



27T 2 d/ 47TK 

1 dlng p 



2tt 2 dl 



0. 



dj/i / 1 / k 2 \\ sin(27TK) . 

- '2-- — +g a h/i + £ "2/1 ( G + G ) + 70/22/3, 



dZ \ 2 \g CT y y 2k 
dy2 ( 1 / k 2 \\ sin(7TK) / ~ 

f - (2 - i + 2/a + ^2/i2/2 + ^2/3 (G + r), 

dG 
~dT 

dG 



= 2(l- 5CT )G + ^(2/ 2 +y 2 ^ 



„ 2^-^)G + .(2/? + 2/ 2 2 ) 



dr 
"dT 

1 dln.g CT _ sin(27TK) / k 2 \ 2 k 2 / ~ 2 jy 2 \ / 2 j/| 



2 (1 - g p ) T + ir {yl + yj) 



4.k K7„- 9 °) yt + 7A G + *)- 9 *{ G + * 

1 rfin. 9p _ _ ffp f| + | + r 2 



2tt 2 dZ ap V 2 2 



A numerical solution of these sets of equations is shown 
in Fig. ([|) and Fig. (||). For these particular plots we used 
2/1,2,3 (0) = 0.01, g a (0) = g p (0) = 1 and G (0) = G (0) = 
r (0) = 0. The RG flows were stopped when any of the 
fugacities reached the value of 1 and the values of the 
other parameters were then plotted at this point. The 
flow equations depend only on the absolute value of k. 

As can be readily checked, the equations always flow 
to strong coupling. Nevertheless, special values of |k| al- 
low us to trace regions with qualitatively different flows. 
Since the RG equations depend only on \k\, those regions 
are mirror reflections on the k = line, the "Toulouse 
line".E£l For k 2 > 3, single spin flip processes are irrel- 
evant (2/1,2 — * OXjust like in the FM single impurity 
Kondo problem.oE 2 ] Moreover, the final flow is indepen- 
dent of the precise value of k , clearly indicating a distinct 
phase of the model. From now on, we will denote this 
phase as region 1. In contrast, spin flips are always rele- 
vant for I At I < 3, but we also encounter a second special 
flow. For |k| = 1, the particle fugacities G and G are 



always the same. There is also a precise balance between 
the "magnetic" (k 2 7 ' g a ) and "electric" (g a ) interactions. 
Consequently, the ground state is a plasma for particles 
of type (m, e, 0), implying that <j) s and 9 S are completely 
disordered. In fact, k = 1 corresponds to the critical 
point of the problem of two weakly coupled LL's. There- 
fore, we can safely identify |k| = 1 as a boundary between 
different phases. For other values of \k\ the interactions 
are screened (g — > or 00) and/or the fugacities have dif- 
ferent flows. It is clear that for 1 < k 2 < 3 (denoted as 
region 2) single-spin-flip fugacities become less and less 
relevant as k 2 — > 3. This suggests a transition region 
from the disordered state at |k| = 1 to the flow of region 
1. We shall call |k| < 1 region 3. In contrast to the 
previous cases, single flips are always strongly relevant 
in this region. The order of relevance of the fugacities 
changes a few times as k is varied in this region. How- 
ever, a particularly simple case occurs in the "Toulouse 
line" (k = 0). 

Even though the renormalization flows are clear and 



(a) Fugacities 



(a) Fugacities 





(b) Interaction constants 



(b) Interaction constants 



Figure 4: RG flow away from half-filling as a function of k. 



Figure 5: RG flow at half-filling as a function of k. 



the special flows were identified, their physical interpre- 
tation is less straightforward. In order to proceed we 
must assign a physical meaning to each particle in the 
gas, from which we can then attempt to determine the 
phase diagram. 



V. EFFECTIVE HAMILTONIANS 

At each RG step we rewrote the problem as a CG. 
Moreover, all the neutrality conditions were preserved 
by the RG step. We therefore can define a quantum 
Hamiltonian that reproduces the CG at each step. This 
effective Hamiltonian allows us to understand the behav- 
ior of the system and, in certain special cases, to infer its 
phase. 

In the dense limit of the KLM, the distance between lo- 
calized spins is of the order of the smallest bosonic wave- 



length available. Therefore, after the first RG step we 
were forced to introduce new entities in the problem. 
Their Hamiltonian form is trivially guessed from their 
definitions, 



Oi 



2 



G 



£ R| r R[ r L] r LJ 



My/ %0.(x) 



S~ (x + S) S~ (x) + h.c. 



G 



pt p pt p 

-RT -RJ. LJ. Lt 



S+ (x + S) S~ (x) + h.c. 



, (18a) 



(18b) 



~ 4r V V F F 

{ r)7 tv}=R,L {ai, a 2 } = hl 

e iy/*Z;* e {x)+4ik F X S z ^ + S z ^ + ^ c .j 18c ) 

where S is a distance of the order of the inverse of the 
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bosonic cut-off (~ a). Both the 0\ and the O2 terms 
involve simultaneous flips of two nearby spins and the 
creation of particles with charges (m, e, c) = (±2, 0, 0) 
and (0, ±2,0), respectively. In contrast, O3 is not re- 
lated to spin flips and generates particles with charges 
(0,0, ±2). It is simple to understand their origin. In the 
original Hamiltonian of Eq. (||) , it is possible to spatially 
resolve the fermion-spin scattering events . As we re- 
duce the bosonic cut-off this is no longer true, and we 
must consider multiple scattering events within the new 
smallest scale, a. There are clear similarities between 
the conduction electron operators in Eqs. ( |i~8| ) and the 
usual backscattering and Umklapp operators. The stan- 
dard picture of the RKKY interaction is that of an ef- 
fective spin-spin interaction mediated by the conduction 
electrons. In light of Eqs. (|l8|), it seems natural to con- 
sider also the opposite point of view: an indirect electron- 
electron interaction mediated by the local spins. The RG 
procedure introduces these composite events in a natural 
fashion. 

The final operator that must be introduced in the ef- 
fective Hamiltonian is a result of the annihilation process. 
Unlike fusion, when a pair is annihilated the zeroth order 
term in an Operator Product Expansion of the bosonic 
fields is a constant. Nevertheless, it is still a function of 
the local spins and must be considered at the last RG 
step in order to establish an effective Hamiltonian. Col- 
lecting all possible pair annihilation terms and expanding 
point-split bosonic operators we get 

O z ~ 4 (g 2 - G 2 ^j S z (x + S) S z (x) 

+ (yl+vl)S- (x + S)S+(x)+h.c. (19) 
It must be stressed that the spin operators in Eqs. flisl ) 



and ( fl9| ) should not be understood as the original local 
spins. Consider the spin history of Fig || as an exam- 
ple. Suppose that the pair flip-antiflip is produced by 
a forward Jj_ term and a backward j\_ term at times r 
and t + 8 within the new renormalization scale. This is 
equivalent to having no flip at all and cannot be distin- 
guished from a particle with fugacity J\. At the operator 
level, this is formally accomplished by summing over all 
possible products of flip operators, expanding the result 
in S ~ a and reordering the Klein factors. The latter 
are actually crucial for the correct final sign (see Ap- 
pendix ^| for details). Thus, we exactly reproduce the 
z-backscattering particle by defining the S z spin at the 
new scale as 



2S Z [x, f) = S+ {x, t + S)S~ (x, t) - S~ (x, t + S)S+ {x, r) . 

A similar calculation can be done for any other possible 
spin history and bosonic operator within a disk of 
radius 6 ~ a. Therefore, the local spins in the effective 
Hamiltonian represent block spins (as in the example 
above) and not the original ones. 
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t t t t_ t 


t t t 
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Figure 6: Two different spin histories at the older RG scale 
that cannot be distinguished at the new one. 



Taking 6 as the lattice spacing at the last RG step 
and collecting all these operators, we find the effective 
Hamiltonian 



Heff = H { 



G - G z 



8rcos 



' ( V 87r 9 P ^c (xj) + Ak F Xj) I S z (xj+x) S z (xj) 



+ 8y 3 sin (yj2-Kg P 4> c (xj) + 2k F xj) sin \yj2-Kg a <j>s (xj)j S z (xj) 

+ 2 yi cos (y/2Trg a (/) s (xjj) + y 2 cos ( y / 27rg P 4> c (xj) + 2k F Xj) e*V K0 ^ x ^ g+ ( Xj ^ 



S + {x j+ x)S (xj) 



GeV f? re ^ ( ^ ) 5+ {x J+1 )S+ ( Xj ) + h. c. 



r 



(20) 



VI. ID ANISOTROPIC KLM PHASE DIAGRAM 



For certain values of K the effective Hamiltonian in 



Eq. (20) is independent of the bosonic fields at the end 



of the RG flow. We will exploit these cases to intuit the 
various phases of the model. 



We start by considering the system away from half- 
filling, where 2/2,3 are irrelevant and T = 0. The RG 
flows are summarized in Table |l|. 

In region 1, the only relevant fugacity is G. There- 
fore, y/8irg a <j) s freezes at it. This reduces Eq. ( |20| ) to the 
anisotropic ferromagnetic Heisenberg model 
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Table III: RG flows for the fugacities away from half-filling. 



(21) 



- G [S* + (x j+i ) S- ( Xj ) + S+ ( Xj ) S- (x j+1 )] 

in its ordered phase (G ~ 1, = 1/2). 

The effective Hamiltonian for the k = line, the 
"Toulouse" line, is also independent of the bosonic field. 
Since the most relevant fugacity is yi, \/2irg a (j) s freezes 
at 7r. This leads to an antiferromagnetic XYZ model in 
an external field 



H, 



eff 



S* {x 3+1 )S z ( Xj )-A yi S x ( Xj ) 




S* (x j+1 )S* ( Xj ) 



s* {x 0+1 )sy ( Xj ). 



(22) 



In this case, the effective spin Hamiltonian exhibits order 
in the XY plane, G ~ G ~ y ~ 1. Nevertheless, this does 
not imply any order of the original spins. As we stated 
before, all our results must be understood in the rotated 
basis of Eq. (|l3|). This ensures that the original model, 
Eq. (||), is still disordered, as emphasized in Ref. [Ts| - 
Therefore, the system is paramagnetic with short range 
antiferromagnetic correlations. Although the "Toulouse 
line" corresponds to a particular case, it seems reasonable 
to extend this assignment to the entire region 3. For one 
thing, because the first term in Eq. (p2|), which drives the 
antiferromagnetic tendency, remains positive throughout 
this region. Besides, the XY disordering terms are the 
dominant interaction in the region. In particular, in the 
|k| = 1 line, the symmetric flow of G and G ensures that 
the z-term vanishes and therefore the order parameter 
(S x ' y ' z ) is still zero. Hence, we propose that the entire 
region 3 is a paramagnetic phase with short range antifer- 
romagnetic correlations. Note that this is not necessarily 
true for other observables, since the flows near |k| = 
and |«| = 1 are qualitatively different. 

There is no simple effective Hamiltonian within region 
2, but the disordering term, proportional to yi, becomes 
progressively less relevant as k 2 — > 3. More importantly, 
the short range z correlations turn from antiferro- to fer- 
romagnetic. Consistent with the identification of region 
1 as a ferromagnetic phase, these two features lead us 



to tentatively identify region 2 as a ferromagnetically or- 
dered phase with unsaturated magnetization of the spins. 

Collecting these results, we conclude that there are at 
least two continuous phase transitions in the anisotropic 
KLM far from half-filling. The first transition, from 
region 1 to region 2 in Fig. ^, reminiscent of the 
Berezinskii-Kosterlitz-Thpuless transition of the single 
impurity Kondo model,c3 separates regions of relevance 
and irrelevance of the single flip process. The effective 
model for region 1, Eq. (|2l|), has ferromagnetic order with 
full saturation of the localized spins. A regime with ferro- 
magnetic order, however, is beyond the present bosoniza- 
tion treatment, since the spin polarization of the conduc- 
tion electrons leads to different Fermi velocities for up 
and down spin electrons. However, the RG flow is still 
able to indicate its existence through the irrelevance of 
single spin flips and the nature of the effective Hamil- 
tonian (pi]). Ferromagnetism is simple to understand in 
the strong coupling limit (\J Z \ Jj_). In this case, ferro- 
magnetic ordering allows the electrons to lower their ki- 
netic energy. In fact, this picture seems ±o survive down, 
to the isotropic case both for J z > OBEl and J z < O.El 
In the FM case, this mechanism is well established and 
is usually called double exchange. However, in the AFM 
chain this simple image of an up-spin electron moving in a 
background of localized down-spins is no longer valid (as 
can be aaal-ytically checked in the Kondo lattice with one 
electrorJl3E3) . In the latter, the objects that lower the ki- 
netic energy are actually the Kondo singlets. Following 
this argument, the total spin per site (electrons+spins) 
would be 5 t z ot = (S z ) — n c /2 = (1 — n c )/2 in the antiferro- 
magnetic case and S z ot = (1 -i-.n«)/2 in the ferromagnetic 
one, as observed numericallyBE 2 ] 

There is another continuous phase transition line from 
region 2 to region 3 in Fig. [?], similar ±o the transition 
of the Ising model in a transverse fieldjEl that separates 
a paramagnetic phase (region 3 of Fig. f?| from a region 
with unsaturated magnetization of the localized spins. 
The magnetization grows continuously up to the border 
of region 1. It is tempting to identify region 2 with simi- 
lar phases with unsaturated moments found in numerical 
studies of both the isotropic FM KLM of DagottO|-ei al.a 
and the isotropic AFM KLM of Tsunetsugu et alH 

The numerical studies of the FM KLM also identified 
a region of phase separation. a We did not find any indi- 
cation of phase separation. We can think of two reasons 
why. First, the coupling constant in that region is of the 
order of the electron bandwidth and therefore bosoniza- 
tion is no longer valid. Moreover, this phase is a com- 
petition between the ferromagnetic tendencies of lower 
band fillings and the antiferromagnetic counterpart at 
half-filling. Since we completely neglect backscattering 
(ultimately responsible for the antiferromagnetism) this 
phase was lost even before we began. 

As we pointed out before, at half-filling we are able 
to include the backscattering terms in the RG scheme. 
We will now consider this case. The first result from 
the RG flow is that O3, whose bosonic part is identical 



12 




H, 



eft 



G — G ) + 4r 



Figure 7: Phase diagram of the ID anisotropic KLM away 
from half-filling. Regions 1 are fully polarized ferromagnets, 
regions 2 are partially polarized ferromagnets and regions 3 
are paramagnetic. See text for details. 



to an Umklapp term, is always relevant, pointing to the 
presence of a charge gap. The spin sector is more subtle 
and we must consider some special cases. 

Region 1 can be simply analyzed. The interaction pa- 
rameters g a and g p go to zero and the most relevant fu- 
gacity is j/3. Therefore, y / 2irg p (f) c freezes at \. The other 
relevant flows are T and G. The effective Hamiltonian re- 
duces to an anisotropic ferromagnetic Heisenberg model 
in a staggered field, 



H eff ~ E(- 4G2 - 8r ) 52 ^+l)^) 
3 

— GS + {x ]+1 )S- (xj)+h.c 
+ (-ir8y 3 S*( Xj ). 



The staggered field induces Neel order, but the system 
has strong ferromagnetic tendencies. As we move away 
from half-filling, the staggered field becomes progres- 
sively irrelevant and the ferromagnetic effective model 
is reobtained. 

The k — point is once again very special. At the 
end of the RG flow, the effective Hamiltonian is also free 
of the bosonic fields and the most relevant fugacities are 
1)2 and y\, They force ^2ng p (f) c and ^/2ng a 4> s to freeze 
at 7T, suppressing the staggered field in the z direction. 
Thus, the "Toulouse point" effective Hamiltonian is 




S z (x j+1 )S* (&) 



S* (x j+i )S* ( Xj ) 



sy (x J+1 )s» ( Xj ) 



- A{y 1 + (-ir V2 )S x {x 3 ). 

As before this does not imply any order of the original 
spins in the XY plane. From this effective model we can 
see that the \n\ = point is characterized by spin and 
charge gaps and no ordering. 

In summary, at half filling we assign two distinct mag- 
netic phases. Regions 1 and 2 have Neel order in the 
z direction. On the other hand, if we assume that the 
"Toulouse line" features can be extended to the entire re- 
gion 3, we can identify this region with a paramagnetic 
phase. The several changes in the relative flows maybe 
a sign of additional phases as a function of k. However, 
the effective Hamiltonian cannot be so easily solved and 
we are unable to make further progress. 

The KLM at half-filling was studied by Shibata et a/.Q 
By looking at the strong coupling limit, they were able 
to find five distinct magnetic phases, which they argue 
survive down to weak coupling: two Neel phases (FM 
and AFM), a planar phase (the triplet state with S z = 
0) , a Haldane phase and a Kondo singlet (paramagnetic) 
phase. Because of the relevance of backscattering and 
the condition J h z -C 1, a direct comparison between the 
RG flows and the available numerical results is restricted 
to k ~ 1. This neighborhood has no simple effective 
Hamiltonian and we are unable to make direct contact 
with the numerical results. We can point out, however, 
that the strong coupling flow of g a is an indication of the 
opening of a spin gap, though this is less certain because 
of the difficulty of analyzing the effective Hamiltonian. 
This possible spin gap is compatible with the Haldane 
type phase at J z < and the Kondo singlet phase at 
J z > obtained in Ref. |lj. As we dope the system away 
from half-filling the backward-scattering terms become 
irrelevant, and a direct comparison with the numerical 
data becomes more feasible. 

As a final illustration of the useful ness of the Coulomb 
gas mapping, we develop in Section VII its application 



to a related yet simplified model of spins and fermions: 
the Ising-Kondo chain. Its simplicity makes it a more 
pedagogical example of the formalism. 



VII. THE ISING-KONDO CHAIN 



The Ising-Kondo model, 
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was proposed by Sikkema ei a/H3 as a model for the weak 
antiferromagnetism of URu2Si2. Here, we will consider 
the one dimensional version of this model and apply the 
same methods that we used in the Kondo chain. Using 
bosonization and disregarding the backscattering terms, 
the Hamiltonian simplifies to 



h = h q + J2 (*) sz (*) - ys x 



(i), (23) 



where the coupling constants were rescaled by the Fermi 
velocity as before. Eq. (|23)-is identical to the co-operative 
Jahn- Teller HamiltonianO The lower symmetry of the 
model allows us to foresee that the sign of J is irrelevant 
to the physics. It is also a well known result from the co- 
operative Jahn- Teller problem that the strong coupling 
limits J » 1 and i/ > 1 show easy axis order in the z 
and x directions, respectively. 

Exactly as in the KLM, we can proceed by going to 
a path integral formulation with bosonic coherent states 
and the local spin S z basis. After tracing the bosonic 
fields and integrating by parts the spin variables, the 
Coulomb gas that follows has only one breed of parti- 
cles (to, 0,0), subjected to the neutrality condition 1 of 
Section II. To mimic our previous notation we define 

g = Assuming y <§; 1 , the RG equation can be de- 

rived in a similar fashion. They correspond to the stan- 
dard Kosterlitz-Thouless equations, 



dy 

dl 
d\ng 

dl 



u 4 = m-9)y, 



-gy 



For g > 1, spin flip processes are irrelevant (see Fig. [|, 
region 1). In the Jahn- Teller language this corresponds 
to a ferrodistortion of the y <C J fixed point. On the 
other hand, for g < 1 spin flips are relevant, y — > oo and 
|, region 2). We can find an effective 
led light on the physics in this regime. 
Indeed, we could have applied the rotation in Eq. (13) to 
the original Hamiltonian to get 



g — > (see Fig. 
Hamiltonian to s 



H = Ho - -e 



S+ (a 



h. c. + s. r. t. 



where "s. r. t." stands for "short range terms". The 



d "displaced" in 
This rotation is 



operators in this rotated basis are ca 
the co-operative Jahn- Teller literature, 
equivalent to the integration by parts of the S z variables 
in the time direction, as we saw. By taking now y — > oo 
and g — > 0, the effective Hamiltonian is simply a magnetic 




Figure 8: Phase diagram of the Ising-Kondo chain. In Re- 
gion 1 the transverse field is irrelevant, while in Region 2 it is 
relevant. See text for details. 



field in the x direction acting to order the local spins. 
Unlike in the KLM, the original spins are also ordered in 
the x direction since 9 freezes at the value of zero. The 
transition is continuous and of the Kosterlitz-Thouless 
type. 



VIII. DISCUSSION AND CONCLUSIONS 

We have proposed in this article what is the natural 
extension to the one-dimensional lattice of the highl^-sufc. 
cessful approach of Anderson, Yuval and Hammar£3Ej'Ell 
to the single impurity Kondo problem. The mapping to a 
Coulomb gas is made specially easy by using bosonization 
methods and particularly subtle developments demon- 
strate the importance of a careful consideration .of Klein 
factors, so often neglected in most treatments.Ea Since 
bosonization relies on the linearization of the conduction 
electron dispersion and is appropriate for the analysis 
of the long-wavelength physics it is never quite obvious 
how far it can be taken in its application to lattice sys- 
tems. However, motivated by its success in the Hubbard, 
Heisenberg and other models, it is reasonable to attempt 
a direct comparison of our treatment to the phases of the 
anisotropic Kondo lattice model. 

One of the hardest tasks in our treatment is the ex- 
traction of physical information from the effective mod- 
els we obtain after several rescaling steps. Some spe- 
cial lines in the phase diagram can be more confidently 
analyzed, but as is common in RG treatments, we are 
then forced to attempt an extrapolation to other regions 
based on continuity arguments. This is specially true in 
our case, where most of the flows are towards strong cou- 
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pling. Given these caveats, however, the overall topology 
of the phase diagram away from half-filling isrfpwipatible 
with the known phases of the isotropic modelBtHl The ex- 
tension of these studies to the anisotropic case would be 
highly desirable. At half-filling, the method itself limits 
its application to the J z <C t region. Unfortunately, this 
is one of the regions where the effective Hamiltonian is 
hard to solve and we are not able to explore the rich phase 
diagram obtained in Ref. |lj. Nevertheless, we do find a 
charge gap at half-filling throughout the phase diagram, 
which seems compatible with the numerical results. The 
question of the spin gap is less clear but our results are 
also compatible with what is known numerically. 

We would also like to try to make contact with previ- 
ous studies of the Kondo lattice model in one dimension 
based on the use of Abelian bosonization. fethe impor- 
tant work of Zachar, Kivelson and EmeryEj where the 
rotation of Eq. ( |l3| ) is first used, the highly anisotropic 
"Toulouse line" (k = 0) is analyzed in detail. One of 
their findings is the presence of a spin gap in the spec- 
trum away from half-filling, which also appears in our 
effective Hamiltonian. At half-filling, they also find spin 
and charge gaps, which seem compatible with our re- 
sults. They also point out that the metal-insulator tran- 
sition-as n — > 1 is of the commensurate-incommensurate 
type.Eil In our treatment, the commensurability condi- 
tion 4fcpa = 27rn = 1 for the relevance of the backward 
scattering terms, the same as in the Hubbard model, is 
a strong indication that the transition is indeed of this 
type. 

Honner and Gulacsi have also investigated the spin dy- 
namics of the isotropic Kondo chain.E3E3 After smearing 
out the discontinuity in the commutation relations of the 
bosonic fields, they replace the latter by their expectation 
value in the non-interacting ground state and write an ef- 
fective Hamiltonian for the localized spins. This Hamil- 
tonian can then be treated numerically and the phase 
diagram determined. This procedure requires the fitting 
of the smearing length scale to numerical results. One 
of the advantages our treatment brings to the problem is 
the ability to do the full analysis analytically and without 
any a priori assumption about the boson dynamics. In 
fact, the Coulomb gas mapping treats spins and bosons 
on the same footing. Besides, no fitting to numerical re- 
sults is necessary. A discrepancy between our results and 
those of Honner and Gulacsi is the partially polarized FM 
phase we find at J z < 0. In their treatment, a PM phase 
is found instead. It would be interesting to extend their 
treatment to the anisotropic case for a fuller comparison. 

Recently, ZacharEl conceived an alternative approach 
to the KLM in the rotated basis. He used a particular 
example of the rotation in Eq. ( |l^ ) 

and treated the KLM in a self-consistent mean-field ap- 
proximation. This approach led him to predict three dif- 
ferent phases in the AFM KLM as well. The first region 
is controlled by the paramagnetic "Toulouse line" fixed 



point. In the rotated basis, this phase is characterized by 
(S z (x)) =0 and (S x ^ 0), precisely as we find in region 
3 of Fig. |. Another phase has (S z ) ji and (S x ) = 0. 
In this case, the system exhibits ferromagnetic order in 
the original basis, and therefore could be identified with 
region 1 of Fig. [?]. Finally, embedded between these two 
phases, he also finds a third intermediate region, which 
he identifies as a "soliton lattice", with (S z ) ^ and 
(S x ) 7^ 0. It is tempting to associate this intermediate 
phase with region 2 of Fig. 0. However, Zachar proposes 
a different description calling region 1 a "staggered liquid 
Luttinger liquid", whereas we find it much more natu- 
ral to associate (S z ) ^ with ferromagnetic order. He 
also conjectures that region 2 does not exist. Finally, 
he argued that all transitions are first order and of the 
commensurate-incommensurate type, while we find them 
to be continuous. 

In conclusion, we have presented a flexible treatment 
of a one-dimensional system of spins and fermions based 
on a mapping to a Coulomb gas, which we treat within 
a renormalization group approach. When applied to the 
Kondo lattice model, the method enables us to identify 
its various phases both at and away from half-filling. 



Acknowledgments 

We thank I. Affleck, A. L. Chernyshev, E. Dagotto, 
M. Gulacsi, N. Hasselmann, S. Kivelson, S. Sachdev, 
J. C. Xavier and O. Zachar, for suggestions and discus- 
sions. E. N., E. M. and G. G. Cabrera acknowledge finan- 
cial support from FAPESP (01/00719-8, 01/07777-3) and 
CNPq (301222/97-5). A. H. C. N. acknowledges partial 
support provided by a CULAR grant under the auspices 
of the US DOE. 



Appendix A: NEUTRALITY CONDITIONS FOR 
THE KLM 



Neutrality conditions are common in Coulomb gas for- 
mulations of quantum problems. In the simplest applica- 
tions, these conditions impose that the overall cjhacge is 
zero as, for instance, in the sine-Gordon model.E30 As 
applied to our case this condition reads 



y^^m{xj,n) = 2J e (xi,Ti) = ^2c(xi,Ti) = 0. (Al) 



They are the mathematical expression of the condition 
for the bosonic correlation functions not to vanish in the 
thermodynamic limit and they also ensure the overall 
cancellation of the Klein factors. However, in the KLM 
the presence of both spins and bosons leads to more strin- 
gent neutrality conditions than in other problems. There- 
fore, besides Eq. (Al), there are two additional restric- 
tions. 
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The first one comes from the impossibility of perform- 
ing two consecutive upward spin flips on a given localized 
spin-i site. Since, from Eq. @ the m = ±1 variable 
gives the direction of a spin flip, it follows that m must 
alternate in time. This condition is also present in the 
Coulomb-gas formulation of the single impurity Kondo 
problemE3 As a direct consequence of the alternation of 
the charge m and the periodic boundary conditions in 
imaginary time, we obtain the first "strong" neutrality 
condition: the total charge to at a given spatial position 
is zero J2i m ( x fixed, T i) — 0. This gives condition 1 of 
Section II, whereas condition 2 is already contained in 
Eq. flO). 



The second additional restriction is slightly less ob- 
vious. From Eq. (Q), we see that each contribution to 
the partition function has a prefactor sign that depends 
on a string of Klein factors and S z operators, the latter 
coming from z backward-scattering events generated by 
H b z of Eq. (^) . The neutrality condition we will derive 
comes from the cancellation of terms with identical ab- 
solute values but with opposite prefactor signs. This will 
finally lead to condition 3 of Section II. We will now 
consider different cases separately. 

Let us first focus on the contributions to the partition 
function coming from terms with forward scattering only 
(Eqs. (|3b"l) and (|3rj)). If there are no spin flips, then the 
prefactor is obviously positive. When there is a pair of 
opposite flips (the only possibility allowed by the neutral- 
ity condition 1), then, because of the overall neutrality 
condition 2, the Klein factors cancel 

^jJcti Fj]a-2 Fr]<7 2 Ft)<ti = Ij 

preserving the positive sign. Consideration of configu- 
rations with additional pairs of flips leads to the same 
cancellation of Klein factors. 

Next, we look at contributions generated by 
ii^(Eq. (|3ej)) only. By considering again increasing num- 
bers of pairs of opposite flips as in the previous paragraph 
we arrive at an analogous cancellation of Klein factors. 

Moving on now to contributions coming from H\ 
(Eq. j3d])), we first consider the possibility of no spin 
flips. In this case, integrating out the bosonic modes, the 
contribution to the partition function is 



a ±2ik F Axi 



S z (t)S z (j), 



(A2) 



where the Klein factors also cancel nicely. Tracing over 
the spin variables leads to no contribution to the parti- 
tion function sum, unless i and j have the same space 
coordinate. What happens for a higher number of inser- 
tions of Hi ? For the general case of N particles coming 
from Hi, the contribution to the partition function will 
be 



LA 



B, 



Z ~XA 



Figure 9: A possib le sp in history and the diagrammatic rep- 



resentation of Eq. (A3). In the figure, A Z B± = B±C Z , where 



A b z = F^F^ain), B A 



F^F^i and CI 



F^Fvatj (r 2 ) 



Note how each S z insertion comes with a corresponding 
c charge. Thus, it is simple to show that tracing over 
S z (i) leads to the condition of having an even number of 
particles of charge c at each spatial coordinate. Moreover, 
the reordering of Klein factors leads to their complete 
cancellation. 

In order to generalize this result to a configuration with 
an arbitrary number of spin flips let us assume initially 
that there are only flips of one kind: either H± or Hj_. 
By using the identity (using Pauli matrices instead of 
spin operators) 



'4"T 2 



cr 2F^ cri F t f, (J2 F^ cr F ucr , 



(A3) 



with v 7^ i] and o\ ^ 02, it is easy show that an Hi 
insertion on one side of a domain wall can be moved to 
the other side with a sign change (see Fig. ^). Now con- 
sider, for example, a pair of particles generated by Hi as 
before. When there is a pair of flips lying along the time 
line, we can move the Klein factors and S z through the 
domain walls with the identity above and cancel them 
out. Therefore, our previous result, obtained without 
the flips, remains valid. This can be generalized for any 
number of flips. Finally, we must consider the possibility 
of having flips coming both from H± and H\ . For this 
we note that a flip from H J ± and a subsequent opposite 
flip from Hj_ "fuse" in a way which is precisely equivalent 
to the insertion of a single Hi particle (Klein factors, S z 
operators and all). Therefore our previous conclusion is 
valid in this case as well: there must be an even num- 
ber of c charges (not necessarily neutral) at each space 
coordinate. 

Finally, we have so far considered insertions along one 
imaginary time line only, which is not the general case. 
Nevertheless, because there is always an even number of 
Klein factors in each time line, we can always reorder 
them so as to group together contributions from indi- 
vidual time lines without introducing additional signs. 
Then, the previous analysis can be used to prove the 
global cancellation of Klein factors and S z operators in 
the general case as well. 

We would like to note that the arguments presented in 
this Appendix indicate a rather surprising precise cancel- 
lation of Klein factors and S z operators, suggesting that 
perhaps there is a deeper underlying symmetry behind 
this result. However, we were not able to find a more 
general symmetry-based demonstration. We also point 
out that, in the problem of a single Kondo impurity in 
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a Luttinger liquid, Lee and Toner0 introduce the same 
kinds of particles denned in the Table |. However, in 
their analysis there is no explicit mention of how to deal 
with the product of Klein factors and the S z operators 
coming from the z backscattering events. We have shown 
that these factors almost miraculously cancel out and do 
not affect the remainder of the analysis of their (or our) 
Coulomb gas. 



Appendix B: ANNIHILATION AND FUSION OF 
PARTICLES 



We now show in detail how the RG procedure leads to 
the annihilation and fusion of charged particles. Consider 
that we initially have a "close pair" with each particle 
having fugacities F\ and F 2 . In the complex notation of 
Eq. (|l0|), the action takes the simple form: 



with: 



/9<T 



zm(rfi) - y/g^e(r)i) 



/9c 



Pij = 



m fa) - V^ e (Vj) + 9 P c fa) c fa) 



-.mfa) + ^fgZefa) 



fg~a~e (rjj) + ——m fa) ) + g p c fa) c fa) 

\j9a 



Suppose the "close pair" particles are at positions I and 
m in space-time. We split the action in three parts: 



+ Pij ln Zi 



Sl = H E + E [aijlnzi 
+ 2ikp c fa) s 

i 

Si = ^2 a a ln zu + pu ln z a 

~\~ ^ ^ C^im ^ %im ~t~ fiim hi Z{ m , 

S3 = aim ln z tm + Pi rn In z lm , 

S\ gives the interaction between the particles which are 
not in the "close pair", S2 the interactions between the 
"close pair" and the other particles, and S3 the interaction 
between the particles belonging to the pair. Finally, we 
define the relative coordinate of the pair as: s = z m — Zi, 
For s ~ di <C 1 we expand the logarithm in s: 



ln Zn 



lnz/i + — . 



Fusion 



If the pair is not neutral, the leading term in the expan- 
sion of S2 is of order zero in \s\. Therefore, we can rewrite 
S2 as giving the interactions between all other particles 
and the new "fused" one. In order to once again write 
the problem in a Coulomb gas form, we must rescale the 
fugacities to accommodate this new particle. Doing the 
integral in S3 



/ 



dse S3 



sin nb 



bh, 



-de, 



where: 



K (e (rim) m fa) + m fan) e fa)) . 



After summing over particle configurations that do not 
contain the fused pair, we get the contribution from fu- 
sion of particles with fugacities F\ and F 2 to the fugacity 
F 3 of this new fused particle 



dF% sin 7r6 f 



dC 



hn 



F1F2 



Annihilation 



If the pair is neutral, the particles annihilate each 
other. In this case, an — —on m and f3u — —fym ■ We 
can expand the partition function contribution in \s\ 



dsdl 



= Si+S 3 



l + \s\B(l,m) + ^-B(l,my + 



where: 



\s\ z a s zu 



The integration over the pair "center of mass" coordinate 
I is 



J dlB (I, m) = 0. 



This is different from the single impurity Kondo prob- 
lem or the dilute limit, where this integral does not van- 
ish. The reason is that, in these cases, the integration is 
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only along the time direction and, therefore, a logarith- 
mic divergence appears. Consequently, the expansion in 
B stops at first order. In contrast, in the dense limit 
the integral is over space and imaginary time, hence re- 
moving this singularity. The first non-vanishing term is 
second otder in B, as in the sine-Gordon and the 2 LL's 
problemEl 



i J dlB{l,m) 



E 



dl 



OtilOtjl 



ZilZjl 



\s\ ZuZjl 
' ZuZjl 



After integration, the first two terms are power law func- 
tions of the distance between the remaining particles of 
the gas. For a sufficiently dilute gas, the most significant 
contribution is given by the last term 



— J dlB (I, to) 2 ~ — 2tt ceuPji In \zij ;| + const 



(Bl) 



It has a simple physical meaning: it gives the "vacuum po- 
larization" coming from the dipole moment of the "close 
pair". The final step in the calculation is to integrate over 
the relative coordinate s 



dse Ss> 



sin nb, 



l m 



blr, 



dL 



(B2) 



where 



km = n{e{rji)m (r? m ) + to (rj m ) e (t?;)) . 

the partition function 



Collecting Eq. (gj and Eq 
contribution after rescaling is 



Appendix C: DETAILED EXAMPLE OF A 
FUSION PROCESS 

In this Appendix we show in more detail how to inter- 
pret the local spins in the effective Hamiltonian of Sec- 
tion [v| after several RG steps. 

Let us focus on the spin histories of Fig. |[ In the new 
RG scale (dashed line), all we know is that the spins at 
times T\ and t 2 have the same orientation. Each process 
compatible with the histories shown in the figure is an 
independent part of the partition function. For definite- 
ness, let us assume that in this position there is a net 
charge (0, 1, 1). At the new scale there are two indistin- 
guishable possibilities to be considered: either there is a 
single particle produced by a term of H%, or there is a 
"close pair" at r and t + St that was fused. 

The effective Hamiltonian strategy is to reconstruct 
the CG at each RG step. Since there is no spin flip 
between T\ and r 2 and there is a net charge (0, 1, 1), the 
operator that performs this task is 

F U F R]S Z & f ) e l v / ^<M^)+V2^^) , (CI) 

where the overbar denotes an operator at the new scale. 

We want to know how to compare the spins at the 
new scale with the ones at the previous scale. In the 
first history of Fig. ^ this is a trivial question. Before 
rescaling, the process had the same form, so 

S z [x,f) = S z (x,t). 

On the other hand, there are four possible bosonic op- 
erators that can fit into the second history case . Let us 
start with the "close pair" 

Fl ] .F Ll Fl l F m S+ (x, r + St) S~ (x, t) 



i ^/TrFg~4> „ ( x , t + 8 r ) + i J |j K 9 „ ( x , r + 8 r ) 
i\Z27Tg p <fi c (x,T) — iJ — k9 s (x,t) 



Point-splitting the bosonic operators, we obtain 



2n—, — — ^2 a iifiji In \z%j I dl 



To complete the RG step, we must sum the charge con- 
figuration of Sithat did not contain the "close par" 



z = e bl 1 - 2ttF 1 F 2 '- 



sin irb, 



bin 



auPjiln\zij\d£ 



Summing over all possible annihilations of pairs of par- 
ticles and reexponentiating, we get the renormalization 
group equations for the Coulomb interaction strengths 
g a and g p . 



J^2tt g t7 tj> s (x,T)+iy/2irg p tf> c (x,T) 



F U F R1 S {x,t + St)S~ (x,t) e 
Another possible pair is 

F U F m F h F Li s ~ (*> T + St ) s+ (*> T ) 

^ y j2-ng p (f) c {x,T J r8T)-\-i^l ^-k6 s (x,r-\-5r) 
iy/27rg cr (j> s (x,r)—i^J k9 s {sc,t) 

Reordering the Klein Factors and point-splitting again, 
we can rewrite this pair as 



F^F m S- (x,t + St)S+ (x,t) 



J.y/'Z-Kg a tf> s (x,f)+i-j2irg p 4> c (x,r) 
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The other two possibilities give the same contributions 
as these ones. If we now identify 

2S z (x,f) = S + (x,t + St)S- (x,t) 
— S~ (x, t + St) S + (x, t) , 



we reconstruct Eq. CI. Note the importance of the Klein 
factors for this identification to hold. 
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